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Abstract 

A new method is presented to obtain the anomalous dimension in 
the solution of the Barenblatt's equation. The result is the same as 
that in the renormalization group (RG) approach. It gives us insight 
on the perturbative solution of the Barenblatt's equation in the RG 
approach. Based on this discussion, an improvement is made to take 
into account, in more complete way, the nonlinear effect, which is 
included in the Heaviside function in higher orders. This improved 
result is better than that in RG approach. 

I. Introduction 



* Mailing address 



1 



Renormalization and the renormalization group (RG) method are power- 
ful mathematical physics tool in quantum field theory and condensed matter 
physics. ^'^ Recently it has been found that there are also important applica- 
tions in dynamical process, asymptotic analysis and perturbation theory.^"^ 
The fractal and chaos behaviors have been studied with the RG approach^. 
The introduction of the intermediate asymptotic idea is very useful. It can 
help us to search for the connection between the similarity, fractal and anoma- 
lous dimension of the RG theory^. But there are still ambiguities in some 
aspects. For example, in quantum field theory and in dynamical theory, it is 
well known that the RG approach can improve the results of the perturba- 
tion calculation. It can remove the divergence and lead to higher leading 
order behavior from a lower order perturbation calculation. It is also clear 
that the divergence and the renormalization are physical requirements for 
some theories. But some of the divergences occur without apparent physi- 
cal reasons. One question involves the improvement that arise from the RG 
perturbative calculations. In this paper we study this problem by analyzing 
a typical example. 

Barcnblatt's equation has many good characteristics. It has strong physi- 
cal background and appears in many different disciplines.® It has been studied 
by applying the method of similarity analysis. It can be solved in terms of 
the parabolic cylinder functions in the intermediate asymptotic region. Its 
perturbation solution is divergent and has anomalous dimension in its asymp- 
totic behavior. It is renormalizable and can be treated with the renormaliza- 
tion and RG approach^'^. But there are still some ambiguities as mentioned 
above. 

In this paper we study improvements on the perturbation solution of the 
Barenblatt's equation by renormalization and RG treatment and, attempt 
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to expose some general principles through the discussion of this example. 
Based on the above analysis, an improved result is obtained by considering 
the nonlinear effect, which is included in the Heaviside function in higher 
orders, in more complete way. 

II. Anomalous dimension calculation by the RG method 

In this section, we repeat some results that are well known, for consistency 

in notation and discussion style. 

We consider the Barenblatt's equation in one space dimension 

1 e 

dtu (x, t) - -dlu (x, t) = -9 {-dtu {x, t)) dlu (x, t) , (1) 

with symmetrical initial condition. It describes a nonlinear diffusive process 

with non- conservative 'mass'. 

Introduce a change of variables 

^^^^"^'^^ = 7^ """"^ T^T{x,t)^t + P. (2) 
Then eq.(l) reads 

dTu{tT)-^d^u{tT)-^dlu{tT) = ^9 [i-dr + ^d^^T)^ dlui^T). 

(3) 

We use an initial condition of u{^, T) as 

It corresponds to the initial condition for u{x, t) in the variables x and t as 
follows: 
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For a general diffusion equation, tliere is a finite solution under this initial 

condition. In this case, the initial condition is singular when approaches 

zero as a 5-function. mo is a finite constant. But it is not so in the case of 

Barenblatt's equation^'^'^. It must be renormahzed to get a finite solution. 

The mo is no longer a constant when the P changes. It becomes infinite when 

the support of the initial state approaches zero. We will discuss that in the 

end of this section. 

Introducing the expansion 

u{^,T) = uo{^,T)+eu^{^,T) + --:, (6) 

one has a sequence of equations: 

2TdTUo{i, T) - ^d^uoii, T) - 9|Mo(e, T) = 0, (7) 

2TdTU,{i,T)-^dM^,T)-dlu,{i,T) = 9 [{-Ot + ^d^)uo{i,T)^ dluo{^,T), 

(8) 

and so on. 

Furtherly, we assume that the initial condition of the uo{^, T) is: 

Then the initial condition of the T) is: 

Un{U^)^Q n>l. (10) 
The solution of the zero-order approximation is: 

Since one has 

[-^T-r^d^ Uo(^,T) = - ^) Uo(^,T), (12) 
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one obtains 



^ (^"^^ + ^^«)«o(e, T)j =e{i- \e\) . (13) 

This is the Heaviside function of the Barenblatt's equation in its lowest order. 
We can write the first-order equation as: 

2TdTU^{^,T) - ^d^u,{^,T) - a|«i(e,T) =9{l- \e\) 5|wo(e,T). (14) 

From the similarity analysis we know that the solution of the Barenblatt's 
equation can be written in a form of variable separation: 

u(tT)^F{T)G(0. (15) 

It is adequate when we are only interested in the behavior of the solution in 
the intermediate asymptotic region. 
Expanding 

F{T)^Fo{T)+eF,{T) + ---, (16) 



and 

we have 
and 

We now take 



G(0 = Go(0+£Gi(0 + ••••, (17) 
Uo{^,T)^Fo{T)Go{0 (18) 
u,{^, T) = Fo (T) Gi(e) + Fi (T) Go{i). (19) 



i^o(r) = -^ and Go(e)=exp(-|). (20) 

The solution in the first-order approximation is 

u^'\i, T) = Fo (T) Goii) + e (Fo (T) G,{i) + F, (T) G^{i)) . (21) 

It tells us that the approximation solution not can be expressed in the form 
of variable separation in finite order of the perturbative theory. But we will 
see that the RG solution is in this form. Thus, we shall pursue along this line 
and accept a corresponding approximation, so that we can use the variable 
separation solution. 
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Inserting eq.(19) into eq.(14) yields 

2T {f[Go + FoGi) - e (FiG; + FoG[) - {f,gI + F^G'[) ^9 {I- F^G^ 

(22) 

where 



^i^.(T), G: = ^G.(0, and = |, 



F[ = —F,{T), G\ = -G,{i), and g", ^ —G,{i), i = 0,l. 

(23) 



Now, we search for the solution of the Fi{T). Integrate this equation with 
respect to the ^ from — oo to +oo and let 

oo 

/ G.iOdC^a. (24) 

— oo 

After simple integral calculation and the application of the conditions 



hrn ^^1(0 = and lirn G,(0 = 0, (25) 

5— >±oo 5— »±oo 



one obtains 



2aFi{T)TFi,{T) + 2V2^TF[{T) + aFo(T) + v^Fi(T) = -^Fo(T). (26) 

Substituting Fq (T) = mo/V^nT and (T) = -Fq (T) /2r into eq.(26), we 
get 

TF[(T) + ^F,iT) + -^Fo(T) = 0. (27) 
The solution of eq. (27) is 

T 
J 

where 



Fi(T) = -F„(r)/Jln^, (28) 



For obtaining this solution, we have made use of the initial condition 

F„(/2) =0 n > 1. (30) 
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The first-order approximation of F[T) is 



F«(r) = Fo(T) (l-aln^), 



(31) 



wfiere 



a = e(3 



e 



(32) 



It is tlie anomalous dimension and is well known in the solution of the Baren- 
blatt's equation in the intermediate asymptotic region by application of RG 
approach^'^. 

In the following discussion, we keep the G(C) in its zero-order approx- 
imation without considering the higher-order corrections. That is to put 
GniO = 0, n > 1. Then the solution of the first order is 



It has been proved that this approximation does not have any effect on the 
function F'^^^T) when we make the RG treatment. We will discuss the space 
behavior in a separate paper. 

For getting a finite but non-zero asymptotic solution, the 'mass' 



not can remain a constant. It is necessary to approach infinity when the 
support of the initial function approaches zero. To obtain a finite solution of 
Barenblatt's equation, the u {x, 0) must be more singular than the 5-function. 
It is different from the diffusive equation in regular cases. Moreover, the 
values of the non-zero order terms are not small, but instead, very large in 
the intermediate asymptotic region. The renormalization is required. It is 
well known that the renormalization and RG approach not only remove the 
divergence but also improve the perturbation expansion solution to higher 
leading order. 




(33) 



oo 




(34) 



— oo 



Let 



mo = mR{i^)Z{i^,l^) 



(35) 



and 



Z(/x,/2) = l + aln 



(36) 
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From the work of Goldenfeld et al,^'^ we know that the result of the renor- 
mahzation and RG treatment is: 

It is easy to check that the initial condition, the boundary conditions and 
the connection condition are all satisfied when the is taken in non-zero 
value. But we must remember that the renormalization mass will approach 
the infinite when P 0, in order to keep the solution finite. Usually we 
require 

X = finite. (38) 

There is an explanation for the improvement of the RG approach.^ But the 
whole case is not very clear yet. We will discuss some aspects of it in following 
sections. 

III. Anomalous dimension calculation by summing F (T) 
to the infinite order 



Now we discuss the improvement to the perturbative solution of the 
Barenblatt's equation in the RG approach. We do this by using the ap- 
proximation suggested in the previous section. We are especially interested 
in the behavior of the intermediate asymptotic region, in which ^ and t are 
far from the zero point but not in infinite. We shall prove that the result of 
the renormalization and RG approach is equivalent to the infinite sum of the 
series of F(r). 

We assume the solution of eq.(3) in a form of variable separation in any 
order. We have known that it is not a suitable form when the approximation 
order n is finite. But, the higher the approximation order is taken, the better 
the proposed form is. As the first step, we keep the dependence on space in 
its lowest-order. Thus we set 

UniC,T) = UT)uoi^,T). (39) 

At the same time, we also neglect a part of the nonlinear effects. We keep 
the Heaviside function in its lowest order as (1 — |^^|) . Therefore we have 
the form of perturbation expansion of the Barenblatt's equation. 

2TdTUr,{i, T) - e%„(e, T) - a|M„(e, T) = 
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n>l. 



(40) 



We will prove that the RG anomalous dimension is just the result of summing 
up the F{T) series to infinite order in the above approximation. 
Inserting eqs.(13) and (39) into eq.(40) yields 

2Tu^{i,T)-^UT) ^e{l- \e\) fn-i{T)dluo{^,T) n > 1. (41) 

The initial condition 



becomes 



n>l 



n>l. 



(42) 
(43) 



When n — 1, one has 



2T«o(e,T)— /i(T) = ^ (l - \e\) a|«o(e,T). 
Integrating it with respect to ^ from — oo to +oo gives 

-P- 



Using fn{l^) = 0, n > 1, the solution of eq.(45) is 

1 T /T^ 



(44) 



(45) 



(46) 



When n = 2, one has 



d 



2TuoiC,T)-^f2iT) ^0{l- \e\) fiiT)dluoiC,T). 
Doing the same integration as we have done to the eq.(44), we obtain 

1 

mo d 



(47) 



(48) 



Due to the identity 



2T^^/i(T) = J didluo{i,n (49) 



-1 

we have 



VrdT 



^/2(T) = h{T)^h{T). (50) 



Making use of /„(r) = 0, n > 1, we obtain 

f2{T) = \[h{T)f . (51) 
Using the principle of mathematical induction, we now prove that 

UT) = ^^[h{T)r . (52) 
Suppose that eq.(52) is true for m = n — 1, that is to say, 

/n-l(^) = ^^^[/l(^)r^ (53) 

We are going to prove that it is true for m — n. Integrating eq.(41) with 
respect to ^ from — oo to +oo gives 

2r^^/„(T) = /„_i(T) / de9^o(e, T). (54) 



Due to eq.(49), one has 



^fn{T)^fn-l{T)^MT). (55) 



Substituting eq.(53) into it gives 



Because /„(/^) = 0, the solution of eq.(56) is 



MT) = ^[MT)]\ (57) 
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It is true for m = 1 and m = 2. Thus we have proved the proposition for any 
m. 

In the infinite limit, we have 



n=l 



(58) 



Inserting eq.(46) into it yields 




n=l 



(59) 



that is 



ii(e,T) = exp^-£/?ln-J uo {C,T) 



Tl 



-T- 



uo iC, T) , 



(60) 



where 



a — e(3 



e 



(61) 



which is just the anomalous dimension given by the renormalization and RG 
approach. 

IV. An improvement on the calculation of the anoma- 
lous dimension 

In the previous section, we have kept the Heaviside function and the 
function G{^) in their lowest order in our calculation of the F{T) to higher 
orders. It has led a partial loss of the nonlinear effect which is included in 
the Heaviside function. Now we show that the calculation of the anomalous 
dimension can be improved by the infinite summing of the function F(T) 
without the approximation of keeping the Heaviside function in its lowest 
order to take account of a part of the lost nonlinear effect. Now we still keep 
the solution in the form of variable separation as eq.(15), but only let the 
function G{S,) in it lowest order. We shall use a trick that makes an analytic 
treatment possible. 

Take the approximation solution obtained in the previous section 




(62) 
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as an initial solution. After inserting it into the Heaviside function, we solve 
the M„(^, T) and search for its infinite sum. Thus we obtain the first iterative 
result. 

We put the first iterative solution into the Heaviside function to solve the 

equation secondly. 

In this way, wc search for the limit of the iterative solution to infinite. If 
the limit exists, we can obtain the solution. 

An important observation is that the form of the solution remains un- 
changed in the iterative process. 

Uo{^,T). (63) 

The only change is the value of the anomalous dimension a'-^^ It simplifies 
the analysis considerably. 
We rewrite a as a^^^ 

a^sP = ^= ^ = . (64) 
v27re 

The a^^^ — 8^5^^^ represents the results after the Mh iteration. 

We still need to solve the problem by summing up the infinite expansion 

series, in each step of the iterative process. 

Putting cq.(63) into the Heaviside function of the nth-order equation, we 
obtain the nth order and [k + l)th iteration equation 

2T^o(e,T)^/r')(T) = 

where the super-index (A;) represents the result of the /cth iteration and the 
sub-index n represents the nth order in the perturbation expansion. The 
initial condition remains in the form of eq.(43). 
When k — l and n = 1, one has 



-a' 



(0) 



2rx.o(e,T)^/«(r) = e [jydT + [i^-) i.o(e,T)j j a|Ko(e,T), 

(66) 
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where fo'\T) = 1 has been used. Because of the identity 



(67) 



one has 

e 



(68) 



Substituting it into eq.(66) and integrating it with respect to ^ from — oo to 
+00 yields 



that is, 



d 



d\nT 



-27(l + 2a(o))exp (-^-^|^) 



-'\/(l + 2a(o))exp 



l + 2o;W^ 



Using = 0, we can obtain the solution of eq.(70) 



where 



A'\t) = /\ In = -/?(^) In 



7(l + 2a(o)) 



When A; = 1 and n = 2, one has 



2Tuoi^, T)^A"\t) = ^ (2a('^) + 1 - e) fi"\T)dluo{C, T). (73) 



(69) 



(70) 



(71) 



(72) 
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Integrating it with respect to ^ from — oo to +00 yields 

^/i^HT)=/«(T)A/f)(T). 
Making use of f^^\P) = 0, n > 2, the solution of eq.(74) is 

Using the principle of mathematical induction, it is easy to prove that 



(74) 



(75) 



1 



(76) 



Suppose it is true for m = n — 1, that is to say 



1 



n-l 



(n-1)! 

We are going to prove that it is true for m = n. Wc know that 

2TuoiC, T)^fi'\T) = e {2a^'^ + 1 - C^) tfi(T)9|«o(e, T). 
Integrating it with respect to ^ from —00 to +00 gives 



(77) 



(78) 



d 

dT 



d 



(79) 



Substituting eq.(77) into it yields 

dT-''' ^ ' (n-l)! 



d_ 
dT 



(80) 



On account of f^^t^) — 0, one can solve eq.(80) and obtain 



ni 



(81) 



We have proved that it is true for m = 1 and m — 2 . Thus the proposition 
is true for any m. 
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Now we sum the infinite series 



oo n 



n=l 



(82) 



Inserting the expression of fi'\T) into it, we obtain 



That is 



where 



n=l 



n! 



v(l) 



and 



(83) 



(84) 



(85) 



The a^^^ is the anomalous dimension corrected to the first-order of the F{T). 
Now, we have finished the first iterative calculation. 

When we make the /cth iteration, the process is the same as the {k — l)th 
one. The only difference is that we will use /S^^-* in place of the (5'^^~^\ Then 
it is easy to prove that 



Mo (C, T) , 



where 



a 



= eP^'^ and P^"^ 



^(1 + 2q;(^-i)) 



(86) 



(87) 



The a^''^ is the anomalous dimension corrected to the kih order of the F{T). 

If the iterative process is convergent, the limit exists. Let limfc_^oo P^''^ - 
p(°°\ which leads to 



P^ 



oo) 



Then one has 



1 + 2e/3M) 

_a(oo) 



(88) 



(89) 
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where 



a 




(90) 



and 



a 



(oo) 



(91) 



^2^g(l+2a(oo)) 



It is obvious that cq.(91) has finite solution, which is the anomalous 
dimension corrected by considering the nonlinear efi^ect of F{T) to infinite 
order. 

The justification of the iterative method and its initial value is that the 
solution is convergent and satisfies the Barenblatt's equation in our approx- 
imation. 

We can compare the improved result with the numerical one given by 
Barenblatt and the result obtained by RG approach as shown in Fig.l. The 
improved result q;(°°) is better than a^"'. 

V. Conclusion 

We have calculated the anomalous dimension by summing the F{T) to 
infinite order. When we keep the Heaviside function and the space dependent 
function G{$,) in their lowest order, the anomalous dimension is the same 
as the RG result. When the approximation of the Heaviside function in 
its lowest order is removed, the value of the anomalous dimension can be 
improved. It is expressed in eq.(91) and Fig.l. 
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